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Summary. — Transverse momentum dependent (TMD) distribution and fragmen- 
tation functions are described as Fourier transforms of matrix elements containing 
nonlocal combinations of quark and gluon fields. These matrix elements also contain 
a gauge link operator with a process dependent path, of which the process depen- 
dence that can be traced back to the color flow in the process. Expanding into 
irreducible tensors built from the transverse momenta p T , we can define a universal 
set of TMD correlators of definite rank with a well-defined operator structure. 

PACS 12. 38. Aw, 12.38.-t, 13.88.+e . 
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1. Introduction 

Transverse momentum dependent (TMD) distribution functions (PDF) and fragmen- 
tation functions (PFF), also simply referred to as TMDs, basically are forward matrix 
elements of parton fields, for instance for quarks 



(1) 



/ -0^z ip< (p\i> 3 {Q)m)\p), 



including a summation over color indices. For a single incoming fermion one would just 
have $ oc (J) + m) . In a diagrammatic approach in principle all sorts of correlators are 
needed, among them quark-quark-gluon correlators defined 



(2) $£ M (p-pi,pib) = 



(2tt) 8 



,i(p-piK e iPi-V 



:p|^(o)a"(»/)^(o|p>. 



(*) Talk presented at Third Workshop on the QCD Structure of the Nucleon (QCD N'12), 
Bilbao, Spain, 22-26 October 2012 
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2 P.J. MULDERS ETC. 

The basic idea is to factorize these hadronic (soft) parts in a full diagrammatic approach 
and parametrize them in terms of PDFs. This will not work for the unintcgratcd correla- 
tors above. At high energies, however, the hard scale provides for each hadron light-like 
vectors P and n such that P-n = 1. For instance n = P' /P-P', where P' is the mo- 
mentum of another hadron obeying P-P' ex s. Using the light-like vectors, one writes 
down a Sudakov expansion of the parton momenta, p = xP + p T + (p-P ~ xM 2 )n with 
x = p + = p-n. In any contraction with vectors not appearing in the correlator, the 
component xP contributes at order -Js, the transverse component at order M and the 
remaining component contributes at order M 2 / 'y/s. This allows consecutive integration 
of the components. Starting from the fully unintegrated result in Eq. 1 one obtains the 
TMD light-front (LF) correlator 



d^-Pd 2 ^ 
(2^ 



(3) ^(x,p T ;n) = J \ n _J T e*** (P\^ 3 (0)MO\P) 



£-n=0 



and the collinear light-cone (LC) correlator 



J ^l JVC m . {0)m )\P) 



or the local matrix element 

(5) *a={P® J (o)MO\F) 



£-n=£ T =0 Or £ 2 =0 
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The importance of integrating at least the light-cone (minus) component p~ = p-P is 
that the nonlocality is at equal (light-cone) time, i.e. time-ordering is not relevant any- 
more for TMD or collinear PDFs [1]. For local matrix elements one can calculate the 
anomalous dimensions, which show up as the Mellin moments of the splitting functions 
that govern the scaling behavior of the collinear correlator $(x). We note that the 
collinear correlator is not simply an integrated TMD. The dependence on upper limit 
<&{x;Q 2 ) = j d 2 p T 3>(x,p T ) is governed by the anomalous dimensions (splitting func- 
tions). One has a a s jp\ behavior of TMDs that is calculable using collinear TMDs and 
which matches to the intrinsic non perturbative p T -behavior [2] . We note that in oper- 
ator product expansion language, the collinear correlators involve operators of definite 
twist, while TMD correlators involve operators of various twist. 

In order to determine the importance of a particular correlator in a hard process, 
one can do a dimensional analysis to find out when they contribute in an expansion in 
the inverse hard scale. Dominant are the ones with lowest canonical dimension obtained 
by maximizing contractions with n, for instance for quark or gluon fields the minimal 
canonical dimensions dim[tp(0)fitfi(^)] = dim[F na (0) F n P(£)] — 2, while an example for 
a multi-parton combination gives dimf^O)^ A"(?y) -0(£)] = 3. Equivalently, one can 
maximize the number of P's in the paramctrization of $y . Of course one immediately 
sees that any number of collinear n-A(rj) = A n (rj) fields doesn't matter. Furthermore one 
must take care of color gaugc-invariance, for instance when dealing with the gluon fields 
and one must include derivatives in color gauge-invariant combinations. With dimension 
zero there is iD n = id n + gA n and with dimension one there is iD" = id" + gA". The 
color gauge-invariant expressions for quark and gluon distribution functions actually 
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include gauge link operators, 

(6) U m =Vexp(-iJ dC M i4"(C)J 

connecting the nonlocal fields, 



(2tt) 3 



(7)$^.( 2; ,p T ;n)= / *' ^ e^< (P\^(0)U m MO\P) 



LF 



(8) *f V\ ^(x,p T ) = j ^'^Y^ eiP<Tr < P ' 5 I F ^^ U m F™{£) U{ m \P,S) 



LF 



For transverse separations, the gauge links involve gauge links running along the minus 
direction to ±00 (dimcnsionally preferred), which are closed with one or more trans- 
verse pieces at light-cone infinity [3, 4]. The two simplest possibilities are U^ = 
UT^ ±oo , UK £i f/S.^ t-i , leading to gauge link dependent quark TMDs $g ' (x, p T ). Which 
correlator is relevant in which process is a matter of doing the (diagrammatic) calculation 
and resummation [5, 6]. For gluons, the correlator involves color gauge-invariant traces 
of held operators F na , which arc written in the color-triplet representation, requiring the 
inclusion of two gauge links £q"o,£] and UL i. Thus even with the simplest links, one has 

already four gluon TMDs $g ' '(x,p T ). 

2. Parametrization of TMDs 

In principle quark and gluon TMDs including a gauge link need to be parametrized 
with a set of PDFs (or in the case of fragmentation PFFs), which just as the correlator 
depend on the gauge link U. For quarks [7, 8] one finds the following set of functions 
depending on x and p\ , 

(9) ^ U \x lPT -n) = |/fW) - flf ] (x,pl) ± r +gW( X , PT )y B 

+ ftgw*^ $ T + htr\ x ,p T ) ^ + iht [u] (x, P i) |}|, 

with the spin vector parametrized as S^ = SlP^+S^+M 2 S L n^ and shorthand notations 
for g[J and /i ls l , 

(10) gf s \x,p T ) = S L g { f2(x,pD - -^-gf^{x,p 2 T ). 

For quarks, these include not only the functions that survive upon p T -intcgration, ff(x) = 
q{x), gl(x) = Aq(x) and h\(x) = Sq(x), which are the well-known collinear spin-spin 
densities (involving quark and nucleon spin) but also momentum-spin densities such as 
the Sivers function fi^{x,p^,) (unpolarized quarks in a transversely polarized nucleon) 
and spin-spin- momentum densities such as giT(x,p^) (longitudinally polarized quarks in 
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a transversely polarized nucleon). The parametrization for gluons, following the naming 
convention in Rcf. [9], reads 

PtSt 

(11) 2xT^W{x,p T ) = -gP ff [U] (x,pl) + g^ 6 -^- f# m (xrf.) 

le 



^* l( * )+ (f-^l) A ""^ ) 



_ e T Pt ,-Lg[U], \ _ t T o T i-e T p T 

2M 2 ^s (X,Pt) m .. 1T 

For TMD correlators time-reversal does not provide constraints as future- and past- 
pointing is interchanged. Depending on the behavior of the Dirac structure, one can 
distinguish T-even and T-odd functions satisfying p™ = ±p u \ where [/* is the time- 
reversed link of U. In the quark-quark correlator / 13 ? and h 1 q are T-odd, in the gluon- 
gluon correlator / 12 ?, h? 1T: h 1 [ and h^ are T-odd. Since time reversal is a good sym- 
metry of QCD, the appearance of T-even or T-odd functions in the parametrization of 
the correlators is linked to specific observables with this same character. In particular 
single spin asymmetries are T-odd observables. 

3. Moments of TMDs 

In order to study the gauge link dependence, it is convenient to construct moments of 
TMDs. The procedure of moments is well-known for the moments of collinear functions. 
For $(x) in Eq. 4 one constructs moments 



^f e** (P|?(0)[^ (iD n ) N i>(0\P) 



LC 



LC 



Integrating over x one finds the connection of the Mellin moments of PDFs with lo- 
cal matrix elements having specific anomalous dimensions, which via an inverse Mellin 
transform define the splitting functions. In the same way one can consider transverse 
moments [4] starting with the light-front TMD in Eq. 3, 



(13) p a T & ± \x,p T -n) = 

^-^- e*** (P|^(0)[/ [ l ±oo] ^ T , ?T]J i??(±oo)^ OOiC] ^(0|P) 



LF 



Integrating over p T gives the lowest transverse moment. This moment involves twist-3 
(or higher) collinear multi-parton correlators, in particular light-cone quark-quark-gluon 
correlator $^ Q (a; — xi,a;i|a;) starting from Eq. 2, and the similarly defined correlator 
<&£){x — x\, x\\x). The particular combinations that are needed in the moments are 

(14) $%(x) = $%(x)-n(x) 

= / (Ixi®d(x — xi,xi\x) - \ dx\ PV — $p a (x — xi,xi\x), 

(15) $%(x)=ir^ a (x,0\x). 
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The latter is referred to as a gluonic pole or ETQS-matrix element [10, 11, 12, 4]. 
The collinear correlators have trivial gauge links bridging the light-like separation and 
thus they have no gauge link dependence. The functions are T-even or T-odd under 
time reversal. While the collinear function $(.t) is T-even, the functions contribut- 
ing to the first moment involve the T-even correlator $>q(x) and the T-odd correlator 
<fr G (x). For the higher moments one finds that the relevant correlators involve similar 
multi-partonic correlators, e.g. for the double weighting one needs the T-even correlator 
<&qq (x) , the symmetrized T-odd correlator $?g G i (x) and the T-even double gluonic pole 

correlator & GG (x). In terms of these functions we get for the moments ^g 1 '^"" 1 ' = 
(16) $M(a;) = *(*), ^ a (x)=^(x) + C^^ G (x), 

(it) ^ ap (x) = *#(*) + 4 a] *?&}(*) +i;cga, c *3s,» > 

c=l 

etc. The gauge link dependence is contained in calculable gluonic pole factors C G , 
etc. For instance the factors C G = ±1. An additional summation is needed if there 
are multiple color arrangements possible for the fields, e.g. the summation over c in the 
double gluonic pole contribution of the quark-quark correlator is needed because one can 
have color structures Tr (GG-0V>) (c = 1) or Tr(tp^)Tr(GG)/N c (c = 2). For the simplest 
gauge links in quark correlators one has C GG 1 = 1 and C GG 2 = 0, but if Wilson loops 
appear in U the latter coefficient is nonzero. 

4. Universal TMDs of definite rank 

The parametrization of quark-quark correlators in Eq. 9 shows that all weightings 
beyond double p T -weighting arc zero. From it we deduce for a nucleon [13] 

(is) & u \x,p T ) = ^ x , P i) + ^i( x y T ) + ^% {x y T ) 



d™ 



(^*k(*,J*) + ^* { W*^>) 



, V^ n [U] Pry -ij , 2 s 
+ Z^ U GG,c ^p ^GG,A X 'Pt)> 
c=l 

where p 1 } is the second rank traceless tensor p 1 } — p l T p J T — \ p\ g? . We refer to the 
correlators in Eq. 18 as TMD correlators of definite rank. Each of these definite rank 
correlators is parametrized in terms of universal TMD PDFs. Using the T-behavior of 
the functions and their tensorial structure the identification is straightforward, with for 
instance fyj, and hj~ appearing in $g- Furthermore, one obtains 

(19) fi(x,pl) = fi(x,pl), gf^{x,p 2 T ) =g 1T (x,pl), ..., 

(20) ft^ ] {x,pl) = C [ ^ft T {x,pl), ht lU] (x,pl) = C G J] h±{x,pl), 

(21) h lT [x,p T )=h 1T (x,p T ) + C GG1 h 1T [x,p T ) + C GG 2 h 1T (x,p T ). 
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Functions like /i and gyr are universal. The T-odd rank-1 TMD PDFs like /ij~ also are 
universal but appear in cross sections with calculable process dependent gluonic pole 
factors. At rank-2, however, one has three universal T-even pretzelocity functions hj~f 
appearing in different, but calculable, combinations in hard cross sections. For a spin 
1/2 target the rank-2 quark-quark correlator $>{qg] — 0> but such a rank-2 correlator 
is relevant for a spin one hadron [13] or for gluon-gluon TMD correlators [14]. The 
analogous treatment for fragmentation functions is simpler. In that case there is no 
process dependence [15]. Phrased in terms of operators, the gluonic pole matrix elements 
vanish in that case [16, 17, 18]. Nevertheless, there exist T-odd fragmentation functions, 
but their QCD operator structure is T-even, similar to the structure of $g. 

The universal p^-dependent TMD functions and correlators show up in azimuthal 
asymmetries, with the azimuthal dependence contained in the tensors in Eq. 18. For 
a given rank m they just are \p T \ m e ±tmip . For the study of the ^-dependence, it is 
convenient to use Bessel transforms [19] of the (— p^,/2M 2 )-moments, 

/ _ 2 \ m /2 /-co 

(22) / (m/2) (x,^) = ( _||j f...(x,p 2 T ) = J db V\^\bJ m (\p T \b)f (m/2) (x,b). 

For quarks one gets multiple color possibilities at rank two, leading to three universal 
pretzelocities. For gluons [14] there are more situations in which multiple color configu- 
rations show up, among them in the study of the gluon Boer- Mulders function h 1 9 (x,p T ) 
corresponding to linear gluon polarization in an unpolarizcd nucleon. 
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